Swarms of particles settling under gravity in a viscous fluid 
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We investigate swarms made of a small number of particles settling under gravity in a viscous 
fluid. The particles do not touch each other and can move relative to each other. The dynamics is 
analyzed in the point-particle approximation. A family of swarms is found with periodic oscillations 
of all the settling particles. In the presence of an additional particle above the swarm, the trajectories 
are horizontally repelled from the symmetry axis, and flattened vertically. The results are used to 
explain how a spherical cloud, made of a large number of particles distributed at random, evolves 
and destabilizes. 



I. INTRODUCTION 

A large number of close but non-touching particles, 
randomly distributed in a spherical volume of a viscous 
fluid and settling under gravity in the same fluid at a 
low-Reynolds-number, forms a swarm, which evolves in 
a remarkable way. The particles circulate and stay to- 
gether for a long time, with the average particle and fluid 
velocity fields the same as in case of a more dense fluid 
dropir— . An initially spherical swarm slowly flattens and 
expands, forming a torus and occasionally leaving single 
particles behind (which, as slower, form a thin long 'tail' 
above the swarm). Then, suddenly, the swarm breaks 
into two (or sometimes more) fragments which repeat the 
same evolution pattern. Dynamics of such a system has 
been extensively investigated experimentally and numer- 
ically, with the use of different methods^ - — . However, 
it is still not clear what is the reason for the observed 
evolution pattern. 

The goal of this work is to construct and analyze sim- 
ple models, which help to understand basic features of 
a swarm evolution. Sec. |TT] contains the guiding ques- 
tions and presentation of the theoretical framework. In 
Sec. IIII1 a new family of periodic solutions is found and 
analyzed, as a model of the swarm. Sec. IIVI illustrates 
how these periodic motions are influenced by the pres- 
ence of an additional particle (to model the tail influence 
on the swarm dynamics). Discussion of the results and 
conclusions are presented in Sec. [V] 



II. PROBLEM 

We aim towards understanding the evolution pattern 
of a swarm of heavy particles settling in a viscous fluid 
at the Reynolds number much smaller than unity. We 
are going to investigate why a swarm slowly deforms (ex- 
pands, flattens, forms a torus) and then suddenly breaks, 
by relating these processes to the circulation of particles 
inside the swarm, and to the occasional loss of the par- 
ticles into a tail above the swarm. To this goal, we con- 
struct and analyze models which contain the characteris- 
tic features of the swarm dynamics (particle circulations, 
their interaction with a tail, deformation of the swarm 



shape and its destabilization), and find out how these 
properties are interrelated. 

The theoretical background is the following. Assume 
that K point-particles, located at move in a fluid 

of viscosity r/ under identical gravitational forces F. The 
fluid velocity v and pressure p satisfy the Stokes equa- 
tions, 

K 

r?V 2 v(r) - Vp(r) = -F^2d(r - n), V ■ v(r) = 0. 



In the reference frame moving with a single particle, 
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where the mobility fi ik is given by the Oseen tensor— 
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-(I + fifcfife) for i ^ k, 
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with r ik = (r t - r k )/r ik and r ik = \rt - r k \. 

III. FAMILY OF PERIODIC SOLUTIONS 

Hocking analyzed oscillations of four point-particles 
settling under gravity in a vertical planes*. His initial 
configurations can be modified to start from point par- 
ticles located at vertices of a rectangle with vertical and 
horizontal sides, all located in the same vertical plane. 

A straightforward 3D generalization results in initial 
configurations of 27V point-particles located at vertices 
of a regular right prism, which consists of 'twin' horizon- 
tal regular A^-polygons, with one particle exactly above 
another one (an example is shown in Fig.QJ. Such config- 
urations lead to periodic relative motion of the particles 
in a vertical plane which includes the center-of-mass of 
the whole group, as schematically indicated in Fig. [1] In 
the reference frame of the center-of-mass, a particle and 
its twin follow the same periodic trajectory. 

The length unit D is twice the initial distance of a 
particle from the symmetry axis, and the time unit is 
8nr]D 2 /F. 
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FIG. 1: Initial configuration of 16 particles (dots). A. Side 
view. B. Top view. Solid lines: planes of periodic motion. 



As an example, we consider evolution of the initial con- 
figuration of 16 particles, shown in Fig. [TJ with c = 0.8. 
The frame of reference is chosen in such a way that the 
z-axis is anti-parallel to gravity, and the xz-plane con- 
tains two periodic trajectories of 2 particles each. All the 
particle trajectories, observed in the laboratory frame of 
reference during t = 25, are plotted in Fig. [2] 

There appear two characteristic stages of the evolution. 



In the first one, for t < 3 (five periods), the particles per- 
form periodic oscillations while settling down. Then, the 
cluster destabilizes, and the motion is not regular any 
more. However, at least half of the particles stay close to 
each other until t = 25, and interact hydrodynamically, 
what is seen as curved fluctuating trajectories. Gradu- 
ally, individual particles separate form the cluster, and as 
slower, are left behind. The second stage of the evolution 
is much longer than the first one. 

In the first stage of the evolution, owing to the rota- 
tional symmetry of the system, all the trajectories are 
planar, and their shapes are identical. 

We now adopt the center-of-mass frame, and investi- 
gate how the shape of the periodic trajectory changes 
with the decrease of c. The trajectory height (just equal 
to c), width w, and its distance to the symmetry axis h/2, 
indicated in Fig. [3l are used as characteristic parameters 
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FIG. 3: Trajectories in the center-of-mass frame for c = 0.1. 

of the trajectories, together with the aspect ratio, 

p = c/w. (3) 

In Fig. [4] we illustrate that the width w and the as- 
pect ratio p of trajectories are getting smaller when c 
decreases down to c=0.2. Their length L decreases, and 




FIG. 2: Trajectories of 16 particles initially located as shown 
in Fig. [TJ with c = 0.8 (side view). 
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FIG. 4: Trajectories in the center-of-mass frame, for large 
values of c = 2.2, 1.4, 0.9, 0.6, 0.2. The smaller c, the shorter 
the trajectory and the smaller its width. 



0.05 



o 




-0.05 - 



FIG. 5: Trajectories in the center-of-mass frame, for small values of c = 0.1,0.05,0.047. The smaller c, the longer and wider 
the trajectory. Note a different scale on each axis. 



a successive trajectory is located inside the previous one. 
When c decreases, the distance h/2 from a trajectory to 
the symmetry axis of the cluster ("radius of the hole") 
increases up to one - the upper limit determined by the 
initial conditions. 



in the log-log scale, see Fig. [5J 




FIG. 6: The trajectory width, w, and twice its distance from 
the symmetry axis, h, versus c, for c > 0.04481. Inset: for 
c > 0.8, the ratio h/w increases linearly with 1/p. 

For small values c< 0.2, trajectories change their shape 
significantly. The length L and width w of trajectories 
increase rapidly when c is decreased by a tiny fraction, 
while h remains close to unity, see Figs. [5] and HO 

As shown in Fig. [JJ the aspect ratio p of a trajectory 
(what can be interpreted as the aspect ratio of a cluster) 
decreases when c is decreased, and for c > 0.3, the slope 
is smaller than for c < 0.2. The inset indicates that the 
trajectory width w becomes very large close to c — c . 

Indeed, we observe that the length L and width w of 
a center-of-mass trajectory, and also period T of the os- 
cillations increase when c decreases. Is there a critical 
value c = Co, where they all become infinite? To check, 
we assume a power law divergence, 

T ~ A/(c - c ) a and to ~ B/(c - c f, for c -> c , (4) 

and search for CQ,a,j3 by plotting T and to versus (c — Co) 




FIG. 7: Aspect ratio p of the cluster as a function of c 
(symbols). Solid straight lines: for larger c (red online), 
p = 0.44c+0.09 and for small c (blue online), p — c— Co. 
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FIG. 8: The period T and the cluster width w versus c — Co. 
Solid lines: T = 0.0013/(c-c ) 1 ' 5 and iu = e /(c-e ). 

In this way, we find the critical value of c, 

c = 0.044788..., (5) 
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and the critical exponents, 

a 1.49 -1.50, ^ w 1.00. (6) 

For c < Co, there is no periodic motions; the swarm 
form the beginning splits into N pairs of particles. The 
comparison of the trajectories without and with periodic 
oscillations is performed in Fig. [9] For c < cq, the slope of 
a trajectory monotonically decreases to a constant value, 
which is smaller than the minimal slope reached at peri- 
odic trajectories. 
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FIG. 9: Two types of the particle trajectories: periodic os- 
cillations (c = 0.05, 0.049, 0.045) and separation into pairs 
without oscillations (c = 0.04, 0.03, 0.02). 



IV. A DROP AND A TAIL 

We will now illustrate how periodic motions, investi- 
gated in the previous section, arc modified in the presence 
of an additional particle. Such a system is supposed to 
model a suspension drop and a tail of particles gradu- 
ally lost from it. It is known that the particles, which 
separate out from a suspension drop, circulate along the 
exterior trajectories up to the top of the drop, and escape 
along the symmetry axisi 

In our model, the drop consists of 16 particles, and the 
tail just from a singlet. The initial configuration of 16 
particles, located at vertices of a regular prism, is shown 
in Fig.[TJ and we take c = 0.9. At t = 0, the 17th particle 
is placed at the symmetry axis of the prism, at a small 
distance zo = 1.25 above its center-of-mass. The question 
is how the trajectories of 16 particles, evaluated in their 
center-of-mass frame (and shown as the middle curves in 
Fig. IU, are modified by the 17^ particle. Can a small 
tail change motion of a large drop? 

The trajectories of the particles from our "drop" are 
shown in Fig. 1101 The particles move in the plane yz 
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FIG. 10: Trajectories in the center-of-mass frame of 16 par- 
ticles, which initially form a regular prism (c = 0.9), with an 
additional single particle on its symmetry axis as indicated. 

of the plot. In the description of the axes, we indicated 
that the reference frame is moving with the center-of- 
mass of the 16 particles which form the "drop". It is 
clear that with time, the shape of the trajectory changes. 
It becomes less high, wider, and its distance from the 
symmetry axis increases. The same features are observed 
when the long-time part of the trajectory from Fig. 1101 is 
compared to the corresponding trajectory in the absence 
of the 17*h particle, see the curve with c = 0.9 in Fig. 2] 
The explanation is that the tail interacts hydrodynam- 
ically with the drop. The additional velocity of a drop 
particle, gained owing to its interaction with the tail (the 
Oseen velocity generated by the single point- force), is 
schematically indicated by arrows in Fig. [TTJ As the re- 
sult, the particles inside the drop are repelled horizontally 
from the drop center, and attracted vertically to the cen- 
tral plane. The closer the tail, the larger the effect. 




; \ 

FIG. 11: Schematic explanation of the horizontal expan- 
sion and vertical flattening of the drop trajectories. Arrows 
(sketch): this part of the drop-particle velocity, which comes 
from its interaction with the tail. 

Similar shape deformation of the trajectories is ob- 
served for spheres rather than points. The results for 
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a drop and a tail made of 8+1 spherical particles are 
discussed in Appendix |A"1 



V. DISCUSSION AND CONCLUSIONS 

A new class of unstable periodic relative motions has 
been found: initially, 2N point-particles form a regular 
prism (two mirror horizontal regular polygons) with a 
vertical-to-horizontal aspect ratio c larger than a criti- 
cal value Co. (For c < Co, the system from the beginning 
separates into pairs of the mirror particles.) Periodic mo- 
tions have been observed for different values of N, with 
a smaller cq for a larger N; for example, cq = 0.017 for 
2N = 32. For larger values of N, it is difficult to inves- 
tigate periodic motions, because they destabilize faster. 
For example, the swarms with w ~ c and N = 16, 32, 64, 
start to move non-periodically after |, |, | of the corre- 
sponding period (T = 1.31,0.52,0.22), respectively. 

The mean settling velocity V s of a swarm with w ~ c, 
averaged over its period, is close to the mean velocity 
Vd = 2NF/(5iri]w) of a fluid drop of the same weight 
and size, with the ratio V s /V d = 0.81,0.88,0.97 for the 
number of particles in the swarm N = 16, 32, 64, respec- 
tively. 

It is supposed that the particle trajectories in the sus- 
pension drop are close to 'a periodic solution'. Particles 
are left behind the drop (as 'a tail') if sufficiently sep- 
arated from these periodic orbits. Such a mechanism 
was confirmed in case of a swarm consisting of 3 parti- 
cles (both point-like and spherical)^— . It was shown 
that 3 close particles circulate together as a swarm for 
a long time before destabilizing into a faster pair and 
a slower singlet, and the interaction time is very sensi- 
tive to the initial conditions 16,17 . The observed chaotic 
scattering was associated with the existence of unstable 
periodic motions which for point-particles have not been 
found explicitlyi 6 -. For 3 spherical particles with a very 
limited class of the initial conditions, such periodic tra- 
jectories of the relative motion were founcU^. For other 
initial conditions, the shape of the relative trajectories 
of the swarm particles indeed resembles the shape of the 
periodic ones^. 

It is expected that the pattern observed for 3 particles 
is generic also for many-particle systems. Periodic mo- 
tions 'hidden' inside a suspension drop can correspond to 
periodic motions of a regular array made of a compara- 
ble number of particles, or a smaller number of effective 
particles. The regular swarm analyzed in Sec. Hill is just 
an example. 

In Fig. [5J we have shown that a regular swarm does not 
break immediately after destabilization. On the contrary, 
most of its particles stay close to each other, and perform 
irregular circulations for a very long time. In this second 
stage of its evolution, the swarm clearly 'stays close' to 
the periodic regular solution from the first stage of the 
evolution. 

Owing to statistical fluctuations, particles which are 



too far from periodic orbits, stop circulating and are left 
behind the drop as a tail above the symmetry axis of the 
drop. This mechanism is consistent with a wide range of 
the drop destabilization times, observed in experiments 
and simulations—, just as it is in case of 3 particles only 16 . 

As demonstrated in Sec. IIVI for a simple model, the 
tail repels particles horizontally from the drop center and 
attracts them vertically to the central horizontal plane of 
the drop. This indicates that larger times t of the drop 
evolution correspond to smaller values of c for the model 
periodic solution of a regular swarm, analyzed in Sec. Hill 
There is an analogy between the dynamics of the drop 
particles at an increasing time t and the dynamics of the 
regular swarm with the decreasing value of the parameter 
c. This analogy is restricted to a limited range of the c 
values, starting from c » w (e.g. c = 2.06 for N = 16), 
and decreasing its value only a little, to match the drop 
aspect ratio p > 0.5 — 0.6 observed in simulations and 
experiments®. For our model, this requirement excludes 
values of c comparable to the critical value Co, see Fig. [7] 

Destabilization time of a suspension drop can be as- 
sociated with the characteristic destabilization time 
of the corresponding periodic orbits. The characteristic 
time scale r* of the drop deformation caused by the par- 
ticle loss seems to be much larger than the characteristic 
destabilization time of the periodic solution. 

In future publications, stability and other basic fea- 
tures of the periodic solutions presented here will be ana- 
lyzed in details, using the methods applied elsewher e 19 ' 20 , 
and related to the corresponding properties of the point- 
particle simulations for suspension drops. Also, it is in- 
teresting to perform a comparison with periodic orbits for 
a regular system of settling non-spherical particles in a 
similar geometry, which have been recently found in sim- 
ulations and experiments 21 . From both fundamental and 
practical point of view, it is important to confirm what 
is the physical mechanism of the drop destabilization. 
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Appendix A: Spherical particles 

Consider now K identical spheres moving in a viscous 
fluid under gravitational forces F. The fluid velocity v 
and pressure p satisfy the Stokes equations with the stick 
boundary conditions at the sphere surfaces. Dynamics of 
the translational motion of the spheres reads 



K 



k=l 



,N, (Al) 



where T"j(t) are time-dependent positions of the sphere 
centers and the mobility matrices fi ik (which depend on 
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relative positions of all the particles) are evaluated nu- 
merically by the multipole expansio n 22 ' 23 with the use of 
the Hydromultipole numerical code^i 

To illustrate periodic gravitational settling of 2N iden- 
tical spheres, we consider the initial configuration of their 
centers at vertices of a cube, with the side equal 1.5c?, 
where d is the sphere diameter (see the upper panel of 
Fig. IT2"j) . Differently from Fig. 0] we show projections of 
the trajectories on the plane of a cube wall (rather than 
the trajectories themselves). Circulations of the spheres 
in their center-of-mass frame are shown in the bottom 
panel of Fig. [H 

Similarly as in Sec. IIV1 we now investigate "a suspen- 
sion drop with a tail" , and we check how does the shape 
of the trajectories in Fig. [12] (drop) change in the pres- 
ence of an additional single sphere (tail), initially located 
above the cube at its symmetry axis, as shown in the up- 
per panel of Fig. [TH] 
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In the bottom panel of Fig. [T3l we show trajectories of 
the drop sphere centers in their center-of-mass frame af- 
ter a long time, when the single sphere has been already 
left up far away. The comparison with Fig [T2"l shows that, 
owing to the hydrodynamic interaction with the tail, the 
width of the drop has become larger, the height smaller, 
and the hole radius has increased, in agreement with our 
findings in Sec. [IV] for the analogical point-particle sys- 
tem. 

Evolution of the system in the reference frame of the 
drop is shown in the movie linked online to Fig. 1131 It 
looks qualitatively the same as for the analogical system 
of point-particles studied in Sec. [TV] Eight spheres circu- 
late inside the drop, while the single tail sphere separates 
from the drop by a rapidly increasing distance. 
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FIG. 12: Model of a drop. Top: initial configuration of 8 
spheres at vertices of a cube. Bottom: Trajectories of these 
sphere centers in their center-of-mass frame. (Projection on 
a vertical plane of a cube wall.) 



FIG. 13: Model of a drop and a tail. Top: initial configu- 
ration of 8 spheres (drop) and one above it (tail). Bottom: 
Trajectories of the drop sphere centers in their center-of-mass 
frame after a time, during which the tail was left up far away. 
(Enhanced online). 
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